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ABSTRACT 

We study the relationship between two major baryonic components in galaxy clusters, namely the stars in 
galaxies, and the ionized gas in the intracluster medium (ICM), using 94 clusters that span the redshift range 
0-0.6. Accurately measured total and ICM masses from Chandra observations, and stellar masses derived from 
the Wide-field Infrared Survey Explorer and the Two-Micron All-Sky Survey allow us to trace the evolution 
of cluster baryon content in a self-consistent fashion. We find that, within r^oo, the evolution of the ICM 
mass-total mass relation is consistent with the expectation of self-similar model, while there is no evidence for 
redshift evolution in the stellar mass-total mass relation. This suggests that the stellar mass and ICM mass in 
the inner parts of clusters evolve differently. 

Subject headings: galaxies: clusters: general - galaxies: elliptical and lenticular, cD - galaxies: luminosity 
function, mass function - galaxies: clusters: intracluster medium 



1. INTRODUCTION 

The baryon content in groups and clusters of galaxies is 
one of the basic quantities that could be used to test models of 
structure formation. The total baryon mass fraction in massive 
clust ers has been used to estimate cosm ological parameters 
(e.g. JWhite e"tai]|1993l:lAllen et al.ll2004l) . The mass of stars 
in galaxies is a fossil record of the star formation and galaxy 
merger history dConrov et al.l 120071) . The relative proportion 
in mass of hot (diffuse intracluster medium, ICM) and cold 
baryons (e.g., stars in galaxies) could be used to gauge the 
strengths of var ious feedback me chanisms and star formation 
efficiency (e.g jBode et"alll2009l) . 

Although many mea surements of the baryon fraction in lo- 
cal clusters exist (e.g.jLin et alj|2003t iGonzalez et alj|2007l : 
lAndreonll2010t iBalogh et alj|201 lb . such studies using large 
cluster samples at higher redshifts are in short supply. The 
first attempt has been made with the X-ra y selected groups at 
z = 0. 1 - 1 detected in the COSMOS field fciodini et al.ll2009l 
hereafter G09). However, the total and ICM masses are not 
well-measured for individual systems, so G09 could not con- 
strain Micm / M sta i- well , especially at the massive end (see also 
iLeauthaudet al.H2011l) . 

Here we present a measurement of the baryon mass frac- 
tion contained in galactic stars and ICM, as well as the rela- 
tionship between Micm and M stal , using 94 clusters at z < 0.6 
- the largest sample with accurate total and ICM mass mea- 
surements to date for this purpose. For the z > 0. 1 systems, we 
use data from the Wide-field Infrared Survey Explorer (WISE; 
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I Wright et al.ll20ldl) survey to determine M star , while for the 
local clusters, w e rely on data from th e Two-Micron All-Sky 
Survey (2MASS. ISkrutskie et al]l2006l) . To ensure uniformity 
of our measurements and to facilitate comparison across the 
cosmic epochs, all quantities are measured within rsoo, the 
radius within which the mean overdensity is 500 times the 
critical density of the Universe p c at the cluster redshift. 

In this paper we have neglected neutral and molecular gas 
in th e baryon budget, as their amount is believed to be small 
(e.g. JChung et alj|2009l) . Our analysis does not include the 
intracluster stars either (as they are well below the surface 
brightness limit of WISE), although they may contribute a 
non-negligible fractio n to the total cluster stellar mass (e.g., 
IGonzalez et alj I2007[). Throug hout this paper we adopt a 
WMAP-5 ( iKomatsu et all 120091) ACDM cosmological model 
where f2 M = 1 - tt A = 0.26, H = 71/z 7 i kms^Mpc" 1 . 

2. ANALYSIS OVERVIEW 

Our analysis relies on the X-ray measurements to provide 
the cluster center, size, and the mass of the ICM. The cluster 
samples used are described in jj2.1| We discuss the way ar- 
tifacts and stars are rejected from the WISE catalogs in jj2.2| 
Throughout this study we only make use of WISE 3.4 pm 
data. For each cluster, we use the background-subtracted total 
flux to estimate the total cluster l uminosity and stellar m ass. 
The method is described fully in iLin et al. (20031 120041) . In 
jj2.3l we outline our modified procedures for analyzing the 
WISE data. 

2.1. Cluster Samples 

We assemble our intermediate redshi ft (z = 0. 1 - 0.6) clus- 
ter sample from | Maughan et all ( 120081 hereafter M08) and 
IVikhlinin et al.l (120091 hereafter V09). both of which are based 
on Chandra observations. Although the 115 clusters from 
M08 are a heterogeneous sample (selected to be targeted 
ACIS-I observations, at z > 0. 1 , publicly available in Chandra 
archive as of 2006), the large sample, analyzed in a uniform 
fashion, is advantageous, and allows us to probe the baryon 
content in clusters of various dynamical states. The Chandra 
subsample of the ROSAT PSPC 400 deg 2 survey presented in 
V09, on the other hand, has a well-defined selection function, 
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FIG. 1 . — The blue and red contours and points show the stellar and galaxy 
loci, combining data from SDSS and UKIDSS DXS. The black points are the 
WISE objects without a match in the 2MASS point source catalog. Only 7% 
of such WISE sources lie along the stellar locus, which suggests that stellar 
contamination in our WISE galaxy catalogs is < 7%. 

but is smaller in size (41 objects). We have taken from these 
studies the X-ray center, ICM and total masfl The total mass 
Msqo is derived from the Y x -M 5 oq relation dKravtsov et alj 
2Q06), where Yx is the product of core-excised X-ray tempera- 
ture T x andMicM- Using clusters common in the two samples 
we have verified that the two analyses yield measurements 
consistent to < 5%. Requiring that the clusters are fully cov- 
ered in the WISE Preliminary Data Release (PDR) footprint 
(see below), and are not affected by nearby bright stars, our 
final sample consists of 45 clusters, with 29 and 16 clusters 
from M08 and V09, respectively. 

As our local (z < 0.1) benchmark , we combine the cluster 
sample presented in iLin et al.l d2004l hereafter L04) with the 
low-z systems from V09. The resulting 49 clusters have ac- 
curate total and ICM mass measurements from Chandra or 
ROSAT, as well as total stellar mass estimates based on the 
near-IR K s -bmd luminosity from 2MASS. 

2.2. WISE Data 

WISE is a satellite mission to survey the whole sky in four 
infrared bands (3.4, 4.6, 12, and 22 pm; hereafter Wl-4), 
with a limiting magnitud^ of 17 (Vega) in Wl. The re- 
cent PDR has made public reduced images and source cat- 
alogs for 57% of the sky. For each cluster we query the 
NASA/IPAC Infrared Science Archive (IRSA) to download 
WISE PDR source catalog out to 5- 10r2oo, where r2oo is a 
proxy of the cluster virial radius, within which the mean over- 
de nsity is 200p r , and is extrapolated from rsoo by assuming 
an iNavarro et al.l (1 19971 hereafter NFW) density profile with 
concentration c = 5. The large area surrounding the clusters 
is used for estimating the local background galaxy count and 
flux. In practice the background is estimated with an annulus 
with the inner and outer radii set to be four and five times r2oo, 
respectively, although the exact choices vary from cluster to 
cluster, primarily depending on the distribution of bright star 

7 For M08 clusters, the measurements are revised with Chandra CALDB 
version 4.3.1. For V09 clusters, all the measurements presented here use 
cir. 2005 calibration; with the most recent calibration, they change only by 
~3%. 

8 http://wise2.ipac.caltech.edu/docs/release/prelim/expsup/sec6_4a.html 



masks. Our results are not sensitive to the exact choice of the 
annulus radii. We only use objects without any known arti- 
facts in Wl and W2. 

Because of its 6" FWHM point spread function (PSF), 
WISE does not provide useful morphological information to 
allow distinction between point and extended sources. In- 
stead, we use 2MASS data to help select galaxies. Specifi- 
cally, for each cluster we download the 2MASS point source 
catalog_| from IRSA, covering the same area as the WISE 
data. We then regard every WISE object with a 2MASS match 
within 2" as point-like and remove them from our catalog. 

To evaluate the efficiency of this approach in star-galaxy 
separation, we use a 1.5 deg 2 region located in the XMM-LSS 
field t hat has deep optical and near-I R data from CFHTLS 
Wide (lGwvnl201ll) and UKIDSS DXS (Lawrence et al.l2007l) 
surveys. In Fig.[T]we show the distribution of the stellar and 
galaxy loci on the g-i vs J-K diagram, obtai ned by cross 
match ing the SDSS star and galaxy catalogs (I Aihara et al.l 
120111) with DXS data. We overlay the 7163 WISE sources 
without a match in the 2MASS point source catalog, and find 
that only 7% of them lie on the stellar locus. This shows that 
using 2MASS we can remove the great majority of the stars. 
For the remaining stars in our catalogs, we remove their con- 
tribution via statistical background subtraction. 

2.3. Estimating Total Luminosity and Stellar Mass 

To describe the evolution of the cluster galaxy popula- 
tions, as well as for the ^-correction, we use a simple stellar 
po pulation evolution model b ased on an updated version of 
the iBruzual & Charlofl (120031) population synthesis software 
package (hereafter the BC model). The model is a single stel- 
lar population of solar metallicity formed at z = 3 with an ex- 
ponentially decaying star formation history (with an e-folding 
timescale r = 0.1 Gyr), with the Kroupa initial mass function 
(IMF), normalized to fit the apparent characteristic magnitude 
m* in the /^-band for massive X-ray c lusters at z = 0.1 -0.8, 
which are taken from lLin et al.l ( 120061 hereafter L06). Using 
the BC model, we follow L06 to construct the composite Wl 
luminosity functions (LFs) in observer's frame in five redshift 
bins. For each bin, we calculate the mean redshift of the clus- 
ters z m , and adjust the ^-correction and distance modulus of 
individual clusters to best represent how the clusters would 
appear i f they were all at z m - We then fit the composite LFs 
with the ISchechterl (119761) function, and show the best-fit m* 
in Fig.|2](top panel). 

For the lower three redshift bins, the Wl m* evolution 
agrees well with the BC model. For the two higher red- 
shift bins, the observed m* is about 0.3 mag brighter than 
the model. To check if this is resulted from the blending of 
galaxies due to the large PSF of WISE, we also construct 
composite LFs in the SDSS z-band and Wl for a sample of 
~ 900 optically-selected cluster s and groups locate d in SDSS 
stripe 82, using the catalog of iGeach et alj (1201 ll hereafter 
Gil). The m* measured for these clusters is shown in the 
middle panel of Fig. [21 along with the BC model predictions. 
Although the measured m* of the SDSS clusters are slightly 
fainter than the predictions, the general agreement of m* in 
both Wl and z-band with the models indicates that blending 
is not the cause of the m* being brighter than the models for 
our X-ray cluster sample at z > 0.4. The galaxies in our X- 
ray clusters are indeed more luminous on average than those 

9 In addition to the standard flags used for clean photometry and to select 
point-like sources, we require the sources to have / < 15.8. 
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FIG. 2. — Upper panel: Evolution of the apparent characteristic magnitude 
m* o f the cluster galaxy LF in Ks-band (blue/solid squares; from lLin et alj 
2006) and Wl (green/solid circles; our X-ray cluster sample). The dotted and 
solid curves show the predictions in K s and Wl of a single stellar population 
formed at z. = 3 with the Kroupa initial mass function and an exponentially 
decaying star formation history, normalized to fit the ^,-band data. Middle 
panel: m* of the cluster galaxy LF in z-band (magenta s quares) and Wl 
(green stars) for the optical clusters from Geach et al. 1 201 1). The two curves 
are the predictions of the same population synthesis model in z-band and Wl. 
Lower panel: evolution of stellar mass-to-light ratio in Wl. 

in the SDSS optical clusters; for example, the mean differ- 
ence of the magnitudes of the brightest cluster galaxies be- 
tween the two samples is about 1 mag. This is not surprising 
given that our clusters are mo re massive than the Gil sam- 
ple on average (see e.g., L04, H ansen et alj|2009l) . With the 
chosen BC model, our evolving m*{z) + 2 flux threshold (see 
below) corresponds to a constant galaxy stellar mass limit of 
2 x 1O 1O M . For this reason, and for sake of consistency, we 
will use the BC model for the ^-correction and m* evolution, 
but will comment below on the impact of the "brightening" of 
galaxies in z > 0.4 clusters on our results. 

For each cluster, we sum the fluxes from objects brighter 
than the smaller of m*(z)+ 1.5 and Wl= 17.0, and subtract 
fluxes from similar objects located within the background an- 
nulus (scaled by the relative area of the annulus and cluster 
region). Assuming all clusters have the same shape of LF 
(with a faint-end power-law slope a = -0.9), which is con- 
sistent with our data (see also L06), we convert the back- 
ground subtracted flux to L tot using an NFW profile with 
c = 5, defined as the total luminosity from all galaxies located 
within r 5 oo that are more luminous than M*(z) + 2. On aver- 
age the conversion from observed to total luminosity requires 
only 15% of extrapolatiorH For reference, at z = we have 

10 Galaxies fainter than M*(z) + 2 would contribute an additional 12% in 
luminosity assuming a = —0.9. 
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FIG. 3. — Correlation between the total stellar mass and ICM mass, both 
measured within rsoo, for the M08 clusters (filled circles), V09 clusters (open 
points), and z < 0.1 clusters (green/squares). To avoid overcrowding the plot, 
eiTorbars for a typical cluster are shown in the lower right corner. We include 
a 10% uncertainty in stellar mass for the assumption that a single value of 
mass-to-light ratio can describe the bulk of galaxies in a cluster. The choice 
of the IMF is the largest systematic uncertainty, at the 30% level. The line is 
the best-fit to clusters at z, ~ (Eqn.|2j. 

T ,wi = l.OlM /L 0)W i- 

For our z < 0. 1 clusters, we measure the total luminosity 
in the X^-band using 2MASS data (following L04); Again 
M s tar = LtotYicXz), using the /G-band Tk,(z) from the same 
BC mod el. Our stellar m ass estimates are consistent with that 
listed in [Lin et al.l (12003b. and are similar to that of lAndreonl 
d2010b and lGonzalez et al.l (120071) . when the intracluster light 
(ICL) contribution is excluded from the latter study. For a 
few clusters at z < 0.1 we have also used WISE data to esti- 
mate M s tar, and found agreement with 2MASS-based values 
to within 10%. 

3. CORRELATIONS BETWEEN STELLAR MASS, ICM 
MASS, AND TOTAL MASS 

In Fig. [3] we show the M stm -M icm correlation, measured 
within r$oo, for the M08 clusters (filled circles), V09 clusters 
(open points), and z < 0.1 clusters (green squares). The M08 
and V09 samples, although constructed using different selec- 
tion criteria, appear to follow the same trend. Together the 
sample spans a factor of 30 in M\cm> allowing us to determine 
the slope of scaling relations well. 

The correlations between the baryon components and the 
cluster mass are shown in Fig. H The lower sets of points 
are the stellar mass fraction (/ star = M star /M 500), and the upper 
sets are the baryon fraction, fi, = (Micm+ Mstar) /A^500> which 
is dominated by the ICM. 

To investigate any redshift evolution among these correla- 
tions, we fit the data with the following form by x 2 minimiza- 
tion, taking into account the uncertainties in both mass mea- 
surements: 

M s tar aM* (1+zp 



ICM*- 

Mstat-cxM^d+zP 
McMfxM^a+zr. 



(1) 



We find that 71 = -0.31 ± 0.22, 72 = -0.06 ± 0.22, and 73 = 
0.41 ± 0. 14; si = 0.62 ± 0.04, s 2 = 0.71 ± 0.04, and s 3 = 1 . 13 ± 
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for simplicity we take 72 = in Eqn.Q]and find that 
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FIG. 4. — The upper panel shows the baryon mass fraction, while the lower 
panel shows the stellar mass fraction. The solid (open) circles are M08 (V09) 
clusters, and the solid squares are z < 0. 1 c lusters. The horizonta l dashed line 
is the cosmic baryon fraction from WMAP (Komatsu et al. 2009). No attempt 
has been made to include the contribution from the ICL. 



0.03. That is, there is no evolution for the M stal -M 5 oo relation, 
weak suggestion of evolution in M stal — Micm, and strong evi- 
dence of evolution in M1CM-M500, which is manifested as an 
offset between the loci of /b for z < 0.1 (solid squares) and 
z > 0.1 (circles) clusters in Fig. |4] This is mainly due to the 
self-similar evolution (SSE) of the ICM, as noted first by V09 
(see the discussion associated with Fig. 10 therein). Let us 
denote Micm oc M^ft , where k = S3 - 1 = 0.13 based on our 
data. Introducing the nonlinear mass scale, Mnl, satisfying 
c(Mnl)D(z) = 6 C 1 .69, where er(M) is the rms flu ctuation 
of lin ear power spectrum and D is the growth factor dPeeblesI 
119801) . we can expect the ICM mass fraction to be the same 
for systems of the same M 5 oo/M NL at different redshifts under 
the SSE. It then follows that Micm oc M'^/M^. We note that 
MnlCz)~ k is consistent with (1 + z) 73 up to z ~ 0.4 or so. In- 
deed, with total mass scaled by Mnl, the scatter reduces from 
12% to 8%. We have confirmed the consistency between the 
observed M1CM-M500 evolution and the SSE expectation by 
repeating the above analysis using M500 based on the T x -M 5 qq 
relation of V09. 

Using all of the clusters we find that 



M stl 



10 12 /^M C . 



= (1.9±0.1) 
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with a scatter of 31%; Without strong evidence of evolution, 
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0.71±0.04 



(3) 

with a scatter of 3 1 %. 

The stellar mass-to-light ratio T is the most important sys- 
tematic uncertainty in our results. We adopt the Kroupa IMF, 
as it gives T at z ~ that is consistent wit h the SAURON 
measu rements of nearby elliptical galaxies (Cappellari et al. 
120061) . Using the Salpeter (Chabrier) IMF, our M star would 
be 34% higher (24% lower). Using a single value of T 
for all galaxies in a cluster is admittedly too simplistic 
dLeaufhaud et al.l 12.01 11) . although we note the spread of T 
in Wl for galaxies of all types is 20% smaller than in opti- 
cal bands. Furthermore, as galaxies in the M* ± 1 magnitude 
range contribute 65% of the total light, our approach is rea- 
sonable, as long as T is only a weak function of galaxy mass 
(and/or morphology). 

Another systematic uncertainty stems from the relative 
mass calibration of the Yx-M^oq relation at different redshifts, 
which is about 5% between z = and z = 0.5 (V09). By boost- 
ing the total and ICM mass by 5% for all clusters at z > 0.45, 
we find that the exponents 71 and 72 become more negative 
(both change by ~ 0.08), giving a weak hint of changing in 
the stellar mass content. 

If only using clusters at z < 0.4, we would have found that 
72 = -0.81 ±0.47 (c.f. -0.06 ±0.22 when using the whole 
sample). That is, the 17 highest-z clusters have a substan- 
tial leverage on the determination of 72. It is also these clus- 
ters that exhibit "brightening" of galaxies with respect to the 
BC model. Should we have adopted the measured m* based 
on these clusters (instead of using the BC model), which is 
equivalent to reducing L tot and adopting a higher galactic stel- 
lar mass limit, we would have obtained 72 = -0.36 ± 0.24. 
We therefore acknowledge the possibility that our results may 
be driven by the behavior of these most massive clusters at 
z ~ 0.5. Ideally we would use a large, volume-limited sam- 
ple that may represent the average cluster properties better, 
which will be carried out in a future publication with the all- 
sky WISE data. 

Possible systematic effects aside, the apparent lack of red- 
shift evolution in M star -M5oo relation is consistent with the 
findings of G09 and L06. In ana lyses that utilize f h to in- 
fer cosmological parameters (e.g.. lAllen et al.ll2004l) . the lo- 
cally derived M stal — M500 relation is usually assumed to hold at 
higher redshifts. Our result directly validates such an assmp- 
tion. 

Our results suggest that, within r 5 oo, the gas and galaxy con- 
tents evolve in different ways; while gas mass grows accord- 
ing to the SSE fashion, the much larger scatter in M star - M500 
and M star - Micm, as well as the much-less-than-unity slopes in 
these scaling relations, suggest a more stochastic growth his- 
tory, which likely involves ti dal interactions to st rip off mate- 
rial from galaxies (e.g., L04. IConroy et al.ll2007l) . In light of 
this, it would be critical to constrain the evolution of the stel- 
lar mass contained in the ICL. The u pcoming Subaru Hyper- 
Suprime Cam survey dTakadall2010l) will likely provide nec- 
essary data in this regard. It is equally important to examine 
the baryon content evolution beyond f-500, both observation- 
ally and theoretically. 
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